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A STRESS  INTENSITY  FACTOR  CALIBRATION 
FOR  CORNER  FLAWS  AT  AN  OPEN  HOLE 

I . Introduction 

The  principles  of  fracture  mechanics  have  been  used  with 
a high  degree  of  success  to  analyze  fatigue  and  structural 
damage  encountered  during  the  operational  life  of  many  struc- 
tures. In  applying  its  principles  to  aircraft  structures,  the 
user  of  fracture  mechanics  must  address  the  problem  of  cracks 
emanating  from  bolt  and  rivet  holes.  A comprehensive  report 
(Ref  5)  on  aircraft  structural  failures  was  completed  in  1971 
and  cited  bolt  and  rivet  holes  as  the  origin  for  33 of  crit- 
ical failures,  by  far  the  largest  single  location.  Realizing 
that  some  of  today's  aircraft  have  over  a million  bolts  and 
rivets,  the  aircraft  designer  must  understand  the  significance 
of  flaws  emanating  from  fastener  holes.  In  addition,  the  Air 
Force  has  recently  adopted  a new  design  philosophy  for  certain 
fracture  critical  parts  which  incorporates  tolerance  criteria 
based  on  the  existence  of  assumed  initial  cracks. 

Using  fracture  mechanics  to  analyze  and  to  design  for 
flawed  fastener  holes  requires  that  the  stress  intensity  fac- 
tor be  known  for  the  flaw  geometry  of  interest.  The  stress 
intensity  factor  is  the  fundamental  parameter  of  fracture  me- 
chanics that  enables  the  calculation  of  the  critical  flaw  size 
for  a given  structure  and  loading  as  well  as  the  estimation  of 
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the  service  life  of  parts  having  subcritical  flaws.  To  date, 
the  three-dimensional  nature  of  flaws  at  fastener  holes  has 
prevented  a complete  fracture  mechanics  solution  for  these 
common  structural  defects. 

The  purpose  of  this  study  was  to  make  stress  intensity 
factor  determinations  for  the  single  geometry  of  a part- thru, 
corner  flaw  emanating  from  a hole  in  a specimen  loaded  in  uni- 
form, cyclic  tension.  In  addition,  this  would  further  estab- 
lish the  fatigue  crack  growth  test  as  a useful  technique  for 
obtaining  stress  intensity  factors  for  difficult,  three- 
dimensional  problems  where  analytical  solutions  are  non-existent. 
The  geometry  of  the  part-thru  corner  crack  and  pertinent  param- 
eters are  shown  in  Fig.  1.  Stress  intensity  factor  calibrations 
for  only  the  end  points  of  the  flaw,  that  is,  points  A and  C, 
were  found.  In  addition,  specimen  thickness,  T,  and  hole 
radius,  r,  remained  constant  from  test  to  test.  This  resulted 
in  corner  cracks  with  a range  in  a/c  (See  Fig.  1)  of  a little 
more  than  1.0  to  about  1.6.  The  most  common  shape  was  approx- 
imately 1.4.  In  sections  to  follow,  the  procedures  and  results 
of  the  fatigue  crack  growth  test  are  described  and  examined. 
First,  a brief  summary  of  background  information  and  previous 
studies  is  presented. 
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Fig.  1.  Geometry  of  a Part-Thru  Corner  Crack 
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II.  Review  of  Information 


Background 

By  applying  the  theories  of  elasticity  to  a flawed  speci- 
men loaded  perpendicular  to  the  crack  surfaces,  one  can  show 
that  the  stresses  local  to  the  crack  tip  can  he  characterized 
hy  a single  parameter,  K,  called  the  stress  intensity  factor. 

K,  in  general,  assumes  the  form 

K = a V^a  MF  (1) 

where  o is  the  far  field  stress,  a is  crack  length,  and  MF  is 
a magnification  factor  that  incorporates  specimen  and  crack 
geometry.  K can  he  thought  of  physically  as  a parameter  that 
describes  the  redistribution  of  stress  due  to  the  introduction 
of  a crack.  K normally  has  the  units  of  psi-VTn.  which  means 
that  MF  is  dimensionless.  In  fact,  it  is  often  desirable  to 
calibrate  (Ref  20)  K in  terms  of  crack  length,  which  can  be 
done  by  solving  Eq.  1 for  MF,  resulting  in  what  is  sometimes 
called  a dimensionless  stress  intensity  factor.  Calibrations 
for  K which  are  given  in  dimensionless  form  (i.e.  in  terms  of 
MF)  can  be  used  to  obtain  K values  for  specimens  of  different 
dimensions  and  varying  magnitudes  of  remote  loading. 

The  determination  of  K is  important  to  fatigue  crack 
growth  analyses.  There  is  a threshold  value  of  K,  say  Kth, 
below  which  a flawed  specimen  can  be  cycled  without  measurable 
crack  extension.  At  the  other  extreme  is  Kc>  called  fracture 
toughness,  a value  of  K at  which  a flaw  will  propagate  unstably. 
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Therefore,  fatigue  crack  growth  must  occur  in  the  range 

to  K . 
c 

If  a specimen  contains  a flaw  (due  to  metallurgical  his- 
tory, fabrication,  improper  fitting,  etc.)  the  size  of  the  flaw 
may  be  too  small  for  incipient  fracture  under  certain  loading 
conditions.  However,  the  flaw  may  be  large  enough  to  grow 
from  subcritical  to  critical  size  under  cyclic  loading.  If  the 
cyclic  loading  is  constant  amplitude  the  crack  growth  can  be 
predicted  by  Paris'  power  law  (Ref  17) 

da/dn  = C(AK)m  (2) 


where  da/dn  is  crack  growth  rate,  C and  m are  experimentally 

determined  constants,  andAK  is  the  stress  intensity  range 

between  which  the  specimen  is  cycled.  Paris'  relationship 

is  the  foundation  of  the  fatigue  crack  growth  test.  C and  m 

are  material  properties,  like  K , which  are  independent  of 

geometry.  Like  K , the  constants  G and  m may  vary  with  testing 

c 

variables  such  as  temperature,  humidity,  and  environment.  Vary- 
ing thickness  from  specimen  to  specimen  also  could  influence  C 

and  m,  just  as  K may  vary  from  a plane  stress  K for  thin  speci- 

c c 

mens  to  a plane  strain  Kc,  called  KIc>  for  thick  specimens. 


Previous  Studies 

Part-thru  cracks  at  holes  present  an  extremely  difficult 
three-dimensional  elasticity  problem,  and  no  exact  solutions 
exist  to  date.  Analytical  and  numerical  approaches  are  few  and 
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primarily  limited  to  estimates  and  models  of  the  problem.  Cruse 
and  Besuner  (Refs  2,4)  have  combined  a boundary- integral  equation 
method  for  elastic  stress  analysis  with  an  influence  function 
method  for  modeling  key  surface  crack  growth  parameters  to  yield 
a technique  for  calculating  K and  subsequent  fatigue  life  of 
surface  cracks  in  areas  of  concentrated  stresses.  Kobayashi 
(Ref  12)  has  extended  an  alternating  method  for  ellipitical 
flaw  analysis  to  solve  the  corner  flaw  problem  in  a quarter- 
infinite  solid.  Finite  element  methods  are  being  applied  in  frac 
ture  mechanics  areas  with  considerable  success.  Cracked  elements 
have  been  created  specifically  for  fracture  mechanics  use,  but  a 
three-dimensional  capability  has  not  been  shown.  Known  surface 
flaw  solutions  have  been  modified  to  yield  estimates  for  the 
corner  flaw  problem  by  Liu  (Ref  13) , Hsu  and  Liu  (Ref  11),  and 
Hall,  Engstrom,  and  Shah  (Ref  8).  In  another  experimental  method 
McGowan  and  Smith  (Ref  15)  used  stress  freezing  photoelasticity 
and  a numerical  method  to  obtain  K's  for  the  corner  flaw  problem. 
The  fatigue  crack  growth  method  of  calibrating  K has  been  used 
most  recently  by  Grandt  and  Hinnerichs  to  study  the  effect  of 
coldworking  holes  (Ref  6).  They  also  presented  some  preliminary 
results  for  surface  flawed  holes  using  PMMA,  a transparent  poly- 
mer. Their  work  and  techniques  formed  a base  of  information 
from  which  this  study  gr ew. 
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III.  Experimental  Program 


The  Method 

With  da/dn  = C(AK)m  as  the  foundation,  the  fatigue  crack 
growth  analysis  can  be  broken  into  three  steps: 

1.  Determine  C and  m experimentally  by  measuring  da/dn 
for  known  K geometries.  (Baseline  tests) 

2.  Measure  da/dn  experimentally  for  the  flaw  geometry  of 
interest. 

3*  Using  the  basic  Paris  relation,  solve  for  AK. 

For  step  one,  several  geometries  are  available,  to  include  the 
Bowie  problem  (Ref  3)  of  a thru  crack  emanating  from  a hole, 
compact  tension  specimens,  and  3 or  4 point  bend  specimens. 

Test  procedures  and  standards  for  baseline  tests  are  given  in 
ASTM-STP-410  (Ref  19)*  The  rationale  for  steps  two  and  three 
is  that  if  the  K parameter  embodies  flaw  geometry,  the  AK  range 
that  is  plotted  versus  da/dn  in  baseline  tests  applies  not  only 
to  the  flaw  geometries  used  in  the  baseline  tests  but  also  to 
any  other  flaw  geometry  of  interest.  This  method  has  the  ad- 
vantages of  being  simple,  inexpensive  (renders  several  data 
points  per  specimen),  and  permits  the  calibration  of  complex 
flaw  geometries.  Its  prime  disadvantage  is  that  it  requires 
data  to  be  differentiated;  therefore,  care  must  be  taken  to 
make  crack  measurements  accurately  and  to  differentiate  properly. 

Materials 

To  enable  direct  observation  of  crack  growth,  polymethyl - 
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methacrylate  (PMMA) , a transparent  polymer,  was  selected  as 
the  test  material.  It  has  the  advantage  of  being  relatively 
isotropic  and  is  also  quite  brittle,  thus  confining  crack  tip 
plasticity  to  a very  small  zone.  PMMA  shows  fatigue  growth 
that  is  macroscopically  similar  to  that  occurring  in  metals 
(Refs  14,15),  and  it  has  been  substantiated  that  the  variation 
in  range  of  stress  intensity  factor,  AK,  controls  cyclic  crack 
growth  over  at  least  a limited  range. 

The  nine  8"  x 14"  specimens  used  in  these  tests  were  made 
from  a 0.72  inch  bulk  plate  of  PMMA.  After  the  holes  were 
bored  and  carefully  reamed,  each  plate  was  annealed  at  100°- 
104°C  for  24  hours  and  then  slowly  cooled  for  24  hours  to  room 
temperature.  PMMA  is  sensitive  to  time,  temperature,  and  en- 
vironment changes.  Also,  thick  components  may  soften  from  heat 
generated  by  higher  frequency  cyclic  loading.  Testing  temper- 
ature (normally  7 0-75°F)  and  testing  conditions  were  not  con- 
trolled but  did  remain  quite  constant.  Cyclic  frequencies 
were  1 or  2 cps  (cycles  per  second);  Arad,  Mukherjee,  and 
Hertzberg  (Refs  1,10,16)  found  that  no  significant  heating  takes 
place  at  less  than  5 cps.  Pitoniak  (Ref  18)  conducted  several 
tests  using  PMMA  at  a frequency  of  3 cps  with  good  success. 

An  MTS  (10,000  lb  capacity)  electrohydraulic  fatigue  testing 
machine  was  used  for  Test  1.  For  all  other  tests  an  MTS  machine 
with  20,000  lb  capacity  was  used. 

Procedures 

A total  of  nine  tests  were  completed:  seven  tests  on  the 
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corner  flaw  geometry  and  two  baseline  tests  to  determine  any 
frequency  effect  of  1 versus  2 cps.  Using  the  data  from  the 
baseline  tests  and  also  data  from  a previous  work  (Ref  6)  the 
constants  C and  m for  PMMA  were  found • Fatigue  crack  growth 
rate  measurements  for  the  corner  flaw  geometry  were  then  used 
to  calculate  the  range  in  stress  intensity  factor  from  Eq  2. 

Previous  baseline  tests  were  conducted  at  a frequency  of 
1 cps  by  Grandt  and  Hinnerichs  using  compact  tension,  3 and 
4-point  bend  specimens.  In  an  effort  to  shorten  testing  time, 
two  baseline  tests  of  1 versus  2 cps  were  conducted  to  deter- 
mine the  effect  of  using  the  higher  frequency.  The  same  spec- 
imens with  a thru  crack  emanating  from  a 0.?5  inch  hole  were 
used  for  the  baseline  tests,  making  them  Bowie  type  problems. 

A sketch  of  a typical  specimen  after  being  epoxy  bonded 
and  bolted  to  grips  for  mounting  in  the  testing  machine  is 
shown  in  Fig.  2.  A small  crack  was  initiated  by  one  of  two  me- 
thods. The  first  was  to  make  a small  notch  with  a file  and 
then  sharpen  the  crack  with  a razor  blade.  The  second,  and 
more  effective  method,  was  to  make  a starter  notch  with  a jew- 
eller’s saw  and  sharpen  the  crack  with  a small  scalpel  blade. 
The  initial  test  was  conducted  on  the  10,000  lb  MTS  machine, 
which  allowed  the  use  of  universal  joints  with  eyebolts  at 
each  end  of  the  specimen.  Subsequent  tests  were  conducted 
on  the  20,000  lb  machine,  and  space  limitations  allowed  the  use 
of  only  eyebolts  at  each  end  of  the  specimen.  This  setup  re- 
quired extreme  care  in  aligning  the  specimen  to  avoid  the  in- 
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Fie:.  2 . Schematic  of  Experimental  Appa 

(From  R 


troduction  of  any  bending  loads.  After  Tests  2 and  3,  it  was 
found  that  one  universal  joint  could  be  included  at  one  end 
of  the  loading  setup,  making  alignment  of  loading  much  easier. 

A small  polariscope  was  used  on  a few  tests  to  check  symmetry 
of  actual  loading  around  the  hole.  While  symmetric  loading 
was  noted,  no  indication  of  bending  could  be  observed. 

The  specimen  was  then  cycled  at  stresses  of  up  to  ?80  psi 
until  a fatigue  crack  started  to  form.  The  load  was  then  re- 
duced to  yield  an  operating  stress  of  590  psi  for  the  corner 
flaw  until  a uniform  crack  front  was  formed.  Cyclic  loading 
was  sinusoidal,  ranging  from  100  lbs  to  3500  lbs  for  a £P  of 
3400  lbs,  with  maximum  deviations  of  ±20  lbs.  This  range  of 
loading  results  in  an  R-ratio  (R  = amin/°max)  of  about  0.029. 
Hence,  R-ratio  effects  were  assumed  negligible  and  were  not 
included  in  calculations.  Km^n  for  the  cyclic  loading  corre- 
sponds to  the  minimum  load  of  100  lbs,  which  is  assumed  to  be 
zero  when  compared  with  the  maximum  load  of  3500  lbs.  Conse- 
quently, AK  becomes  equal  to  K.  AP  of  3400  lbs  versus  P of 

max 

3500  lbs  was  used  for  all  calculations. 

Crack  growth  was  recorded  by  means  of  a 35  mm  camera  with 
bellows.  Using  the  mirror  setup  as  shown  in  Fig.  2,  pictures 
were  made  directly  through  the  end  perpendicular  to  the  crack 
using  TRI-X  (ASA  400)  black-and-white  film.  Lighting  of  the 
crack,  especially  in  early  growth  stages,  was  particularly 
difficult.  A high- intensity  lamp  aligned  parallel  and  as  close 
to  the  camera  as  possible  was  found  to  provide  the  best  illumi- 
nation. 
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Measurements  of  crack  lengths  on  the  35  mm  negatives  were 
made  using  film  readers  that  magnified  the  image  approximately 
ten  times.  Photographs  were  taken  at  varying  cyclic  intervals 
that  allowed  measurable  crack  growth  of  as  small  as  0.005  inch. 

A scaling  factor  for  the  measurements  was  formed  using  the  hole 
diameter  as  a reference  length.  Some  typical  crack  photographs 
are  shown  in  Figures  3 and  4.  Reflections  of  the  crack  and  both 
surfaces  of  the  plate  were  very  troublesome  when  measuring  small 
cracks,  and  this  reflection  can  also  be  seen  in  the  photos. 

The  Bowie  baseline  tests  were  conducted  using  the  same 
procedures  as  above  with  few  exceptions.  The  initial,  thru 
crack  was  made  with  a jeweller's  saw,  and  then  the  specimen 
was  cycled  until  a complete,  thru-the-thickness , fatigue  crack 
had  formed.  Photographs  were  then  made  under  a sinusoidal, 
cyclic  loading  of  100  to  2300  lbs.  Again,  the  R-ratio  was 
approximately  0.044,  and  its  effect  was  thus  neglected. 

A process  that  was  used  by  Grandt  in  a previous  work  (Ref 
7)  was  used  to  differentiate  present  crack  growth  versus  cycles 
data.  Five  successive  pairs  (a,n)  of  data  points  were  consid- 
ered in  a set.  A least- squares  parabola  was  passed  through 
these  points  and  crack  length  calculated  at  four  regular  inter- 
vals over  the  range  of  the  set.  The  growth  rate  was  computed 
by  a standard  least- squares  formula  which  approximates  the 
derivative  at  the  center  of  five  evenly  spaced  points.  The 
set  was  then  advanced  one  point  and  the  process  repeated  until 
all  the  data  was  exhausted.  Slightly  different  least-squares 
formulas  were  used  for  da/dn  near  the  ends  of  the  range  of  data. 
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. 3*  Photo  of  a Corner  Flaw,  Test  9 (14,000  cycles) 
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. 4.  Photo  of  a Corner  Flaw,  Test  9 (27,000  cycles) 
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IV . Experimental  Results 


Frequency  Effect 

The  fatigue  crack  in  both  tests  initiated  slightly  off 
the  centerline  of  the  thickness  and  then  grew  until  the  crack 
front  was  completely  thru  the  thickness.  Then  the  crack  assum- 
ed a stable,  almost  symmetrical  shape  as  shown  in  Figures  5 and 
6.  An  average  crack  length  was  computed  using  the  average  of 

at’  am’  and  ab’  The  crack  length  at  one  of  the  plate's  sur- 
faces is  at;  am  is  the  longest  crack  length  (usually  at  mid- 
thickness); and  a^  is  the  crack  length  at  the  other  surface. 

For  the  1 cps  test,  a^  was  not  used  in  the  average  because  pho- 
tograph quality  did  not  permit  accurate  measurement  of  a. . 

Using  Bowie's  solution  (Ref  3)  for  a thru-crack  emanating 
from  a hole,  stress  intensity  factors  for  the  various  crack 
lengths  were  calculated.  The  log  of  crack  growth  rate  was 
then  plotted  versus  log  AK,  and  the  results  are  shown  in  Fig- 
ures 7 thru  11.  The  least-squares  straight  line  which  best 
fits  each  set  of  data  is  also  drawn,  and  the  plots  clearly  show 
the  linear  relation  of  log  da/d n versus  logAK  described  by 
Eq  2.  Plots  of  the  separate  tests'  data  are  shown  in  Figures 
7 and  8,  and  a plot  of  the  combined  1 and  2 cps  data  and  its 
corresponding  Paris  equation  is  shown  in  Fig.  9.  Both  the  1 
and  2 cps  data  appear  to  follow  the  same  Paris  equation  fairly 
well,  but  perhaps  more  examination  is  needed  before  any  conclu- 
sions about  frequency  effect  can  be  made.  In  Fig.  10  the  least- 
squares  line  which  best  fits  the  1 cps  data  is  compared  to  the 
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Photo  of  Bowie  Crack,  Test  4 (55.000  cycles) 
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Photo  of  Bowie  Crack,  Test  7 (23,000  cycles) 
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Fig.  11.  Baseline  Data,  (138  points) 
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least-squares  best  fit  of  the  2 cps  data.  The  very  small  effect 
of  testing  at  the  two  frequencies  is  now  more  evident.  When 
using  baseline  data  one  usually  must  set  limits  to  define  the 
range  over  which  the  data  will  give  good  results.  The  limits 
that  were  selected  for  this  study  were  those  used  by  Grandt  and 
Hinnerichs  in  Ref.  6;  specifically,  if  da/dn  were  less  than 
0.15(10)-5  inch/cycle  or  greater  than  0.3(10)”-^  inch/cycle, 
stress  intensity  factor  and  magnification  factor  were  set  equal 
to  zero.  These  limits  were  substituted  into  the  Paris  equations 
of  the  separate  1 and  2 cps  sets  of  data,  the  equations  of  which 
are  given  in  Fig.  10,  and  stress  intensity  factors  were  calcu- 
lated. At  the  upper  limit  the  difference  between  the  two  values 
was  less  than  5^»  and  at  the  lower  limit  the  difference  was 
less  than  1$.  This  was  considered  a very  small  difference  for 
stress  intensities  in  the  range  400-1100  psi-Vin,  and  it  was 
thus  concluded  that  tests  could  be  run  at  1 or  2 cps  with  no 
significant  loss  of  accuracy.  For  evaluating  present  data, 
the  96  Bowie  baseline  data  points  were  combined  with  the  226 
baseline  data  points  of  Ref.  6 to  produce  a baseline  data  base 

of  322  points.  The  constants  yielded  from  the  322  points  were 

-20 

C = 0.22138(10)  and  m = 5-591*  Since  a plot  of  the  322  points 
tends  to  make  individual  data  points  indistinguishable,  138 
points  were  selected  from  the  full  data  base  and  are  shown  in 
Fig.  11. 

Comer  Flaw 

Seven  tests  were  conducted  on  the  corner  flaw  geometry. 

Test  numbers  1 and  2 were  conducted  at  1 cps,  and  Tests  3,  5, 


23 


6,  8,  and  9 were  done  at  2 cps.  Plots  of  crack  length  vs. 
cycles,  crack  shape  (a/c)  vs.  a/T,  magnification  factor  vs. 
a/T,  and  stress  intensity  factor  vs.  c are  given  for  each  test 
in  Figures  24  thru  51  in  the  appendix.  (See  Fig.  1 for  a visual 
explanation  of  the  quantities  a,  c,  and  T.)  Crack  lengths  at 
the  surface,  c,  and  along  the  hole,  a,  are  plotted  along  with 
their  associated  stress  intensity  factors,  KC  and  KA.  Also, 
dimensionless  magnification  factors  for  doth  a and  c crack  lengths 
are  plotted,  with  magnification  factor  at  the  hole,  MFA,  being 
given  by 


MFA  = (KA)/(o  Vrra)  (3) 

Magnification  factor  at  the  surface,  MFC,  is  given  by 

MFC  = (KC)/(o  V^c)  (4) 

The  ratio  a/T  is  a dimensionless  quantity,  T being  thickness 
in  inches.  (Note  for  Test  2s  The  odd  portion  of  the  crack 
growth  curve  was  caused  by  the  lack  of  photographs  for  that 
portion  of  the  growth. ) 

Specimen  dimensions  are  generally  a width  of  8 inches, 
thickness  of  0.72  inch,  and  a hole  diameter  of  0.75  inch. 
Actual  specimen  dimensions  for  each  test  specimen  are  given 
in  Table  II  in  the  appendix.  Tests  covered  periods  of  time 
ranging  from  hours  to  four  days.  It  was  found  by  Grandt  (Ref 
7)  that  interruptions  had  no  apparent  effect  on  test  results. 
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Crack  growth  was  characterized  in  early  stages  by  growth 
to  the  natural  crack  shape.  The  term  natural  refers  to  the 
tendency  of  a starter  crack  in  constant  geometry  specimens  to 
assume  about  the  same  a/c  ratio  when  subjected  to  the  same 
loading.  The  a/c  ratio  in  all  cases  was  greater  than  one. 

Early  growth  to  the  natural  shape  can  be  seen  in  the  first 
portion  of  the  a/c  curve  for  Tests  6 and  8,  Figures  4l  and  45. 
Some  of  the  oscillation  in  early  crack  length  measurements  was 
probably  caused  by  the  combination  of  poor  lighting  and  reflec- 
tion, which  produced  a poor  quality  picture.  In  general,  the 
quality  of  pictures  increased  as  the  crack  grew,  resulting  in 
better  illumination. 

It  was  also  observed  in  both  baseline  and  corner  flaw 
tests  that  the  fatigue  crack  front  at  the  surface  (point  C in 
Fig.  1)  lagged  behind  the  crack  front  within  the  thickness. 

This  is  probably  direct  observation  of  the  transition  from  a 
plane  strain  to  a plane  stress  condition  as  predicted  for  a 
thick  specimen.  Some  specimens  were  cycled  to  fracture,  and 
the  fracture  surfaces  showed  no  shear  lips,  indicating  that 
plane  stress  zones  were  very  small  and  that  plane  strain  frac- 
ture had  definitely  occurred. 
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V.  Discussion  of  Results 


The  works  of  McGowan  and  Smith  (Ref  15)  and  Hall  and  Finger 
(Ref  9)  both  noted  a decrease  in  stress  intensity  factor  at  the 
hole,  KA,  and  an  increase  in  stress  intensity  factor  at  the  sur- 
face, KC , as  a/T  reached  higher  values.  McGowan  and  Smith  con- 
jectured that  this  result  might  be  due  to  the  fact  that,  when  the 
crack  is  deen,  a substantial  part  of  the  load  is  transferred  to 
the  side  of  the  hole  opposite  the  crack  due  to  its  greater  stiff- 
ness. Thus,  the  remaining  ligament  between  the  "a"  crack  front 
and  back  surface  of  the  plate  carries  a reduced  load.  Another 
point  could  also  be  considered.  As  the  crack  approaches  the  back 
surface  of  the  plate,  the  ligament  may  become  part  of  the  plastic 
zone  at  the  crack  front,  causing  load  redistribution. 

The  decrease  in  KA  and  increase  in  KC  was,  in  general,  not 
in  evidence  in  present  test  data.  This  is  probably  because  plas- 
ticity effects  were  kept  to  a small  scale  by  using  small  loads 
in  the  fatigue  test.  Pitoniak  (Ref  18)  found  several  values  for 
PMMA  yield  strength  and  fracture  toughness,  averages  of  which 
are  7000  psi  and  1000  psi-Vln,  respectively.  The  present  tests 
used  a maximum  stress  of  about  600  psi,  less  than  0.1  of  the 
yield  strength.  This  stress  yielded  maximum  stress  intensi- 
ties near  the  crack  front  on  the  order  of  100  psi-Vin'»  again 

0.1  of  Kt  . Therefore,  it  can  be  concluded  that  elasticity  was 
Ic 

on  a small  scale,  and  load  redistribution  did  not  occur  at  the 
higher  a/T  values.  In  general,  it  can  be  supposed  that  obser- 
vation of  the  KA  and  KC  changes  are  material  and  load  dependent. 
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Figures  14  thru  23  in  the  appendix  are  plots  of  dimension- 
less magnification  factor  versus  normalized  a/T.  There  are  five 
sets  of  "a"  and  "c"  magnification  factors  for  five  ranges  of  a/c. 
The  data  points  come  from  all  seven  tests  on  the  corner  flaw 
geometry  and  encompass  86^5  of  the  total  number  of  test  points. 
Results  compare  favorably  with  the  preliminary  tests  of  Grandt 
and  Hinnerichs.  Comparisons  are  given  in  Figures  12  and  13  of 
part  of  the  data  with  the  estimates  of  Granf  et.al. , Hsu  and  Liu, 
and  Hall  and  Finger  (Refs  5,11,9)-  A discussion  of  the  three 
estimates  and  how  they  compare  with  present  data  now  follows. 

The  present  data  points  shown  in  the  comparison  plots  are  mag- 
nification factors  for  the  1.35~1*45  a/c  range.  Accordingly, 
points  used  to  plot  the  estimates  were  calculated  using  nominal 
values  of  a/c  of  1.4,  thickness  of  8 inches,  and  a hole  radius 
of  0.37  inch.  Further  comparison  between  the  estimates  and 
present  data  is  given  in  Table  I.  Present  data  given  in  the 
table  are  from  Test  1 where  pertinent  parameters  agreed  quite 
closely  with  the  nominal  values  used  in  the  estimates  calcu- 
lations . 


Gran,  et  al . 

The  basic  solution  for  an  embedded  elliptical  flaw  was 
modified  to  account  for  hole  and  free  surface  effects  and  also 
included  a correction  factor  for  out  of  shape  ellipses.  The 
final  expression  for  stress  intensity  factor  at  the  hole  (KA)  is 


KA 
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Table  I 

Comparison  Table 


Comparison  of  Hole  Magnification  Factor  Estimates: 


a/T 

Present 

Data 

Gran 

Hsu/ 

Liu 

Hall/ 

Finger 

Modified 

Hall/Finger 

• 3 

1.39 

2 . 60 

2.12 

1.22 

1.40 

.4 

1.31 

2.69 

2.12 

1.17 

1.35 

• 5 

1.25 

2.82 

2.12 

1.09 

1.25 

.6 

1.22 

3.01 

2.12 

1.03 

1.18 

.7 

1.19 

3.33 

2.12 

.98 

1.13 

.8 

1.23 

3-90 

2.12 

.94 

1.08 

fo  Difference  of 

MFA  Estimates  from 

Present  Data: 

• 3 

— 

+87 

+53 

-12 

+1 

.4 

-- 

+105 

+62 

-11 

+3 

.5 

-- 

+126 

+70 

-13 

0 

.6 

-- 

+147 

+74 

-16 

-3 

• 7 

-- 

+180 

+78 

-18 

-5 

.8 

+217 

+72 

-24 

-12 

Comparison  of  Surface  Magnification  Factor  Estimates: 


a/T 

Present 

Data 

Gran 

Hsu/ 

Liu 

Hall/ 

Finger 

Modified 

Hall/Finger 

•3 

1.58 

1.34 

1-57 

1.44 

1.66 

.4 

1.46 

1.21 

1.44 

1.38 

1-59 

• 5 

1.41 

1.13 

1.32 

1.29 

1.49 

.6 

1-35 

1.06 

1.24 

1.22 

1.40 

• 7 

I.36 

1.02 

1.18 

1.16 

1.33 

.8 

1.36 

•99 

1.14 

1.11 

1.27 

% Difference  of  MFC  Estimates  from  Present  Data: 

• 3 

-- 

-15 

-1 

-9 

+5 

.4 

-- 

-17 

-1 

-5 

+9 

• 5 

-- 

-20 

-6 

-9 

+6 

. 6 

-- 

-21 

-8 

-10 

+4 

• 7 

-- 

-25 

-13 

-15 

-2 

.8 

-- 

-27 

-16 

-18 

-7 
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Here  5 is  the  complete  elliptical  integral  of  the  second  kind. 
A plot  of  $ values  versus  a/c  is  given  in  Ref  11.  This  yields 
a magnification  factor  at  the  hole  (MFA)  of 


™ ‘ I *an  § > 


(6) 


For  the  stress  intensity  factor  at  the  surface  (KC), 


KC  = 1.120,  ( 1 + i(-X-)2  + ^(^)4  ) 


2xr+c'  2'r+c 


(?) 


and  magnification  factor  at  the  surface  (MFC)  is  given  by 


MFC  = 


1.12  / 1 . 1/  r s2  . 3/  r » 

( 1 + 2(F^}  + ) 


2 r+c 


(8) 


where  b was  set  equal  to  c. 

It  can  be  seen  from  Table  I that  the  estimate  for  magni- 
fication factor  at  the  hole  is  considerably  larger  than  present 
data  at  a/T  = 0.3  and  diverges  with  increasing  a/T;  therefore, 
this  estimate  is  not  included  in  the  plot  of  Fig.  13.  The 
estimate  for  the  surface  magnification  factor  is  much  closer 
to  present  data,  ranging  from  15#  to  27#  less,  with  an  average 
of  21#  less. 

Hsu  and  Liu 

Stress  intensity  expressions  for  truncated  elliptical 
cracks  at  the  comer  of  a quarter-infinite  solid  were  derived. 
These  were  modified  to  account  for  the  influence  of  the  hole 
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and  the  front  free  surface  to  yield  the  following  expression 
for  magnification  factor 


MF  = B §1 

$ 


(9) 


B is  Bowie's  factor  (Ref  3)  which  accounts  for  the  influence 
of  the  hole  on  a full-thru  crack,  and  Mj/S  is  a factor  that 
accounts  for  the  front  surface  influence  and  also  flaw  shape. 
Plots  of  are  given  in  Ref.  11.  For  MFA,  B = 306  and 

Mj/a  is  also  constant,  equaling  O.63.  Then, 


for  all  values  of  a/T.  For  MFC,  M^/5  is  constant  and  equals 
0.86;  but  B varies  with  distance  from  the  hole.  Therefore, 


Table  I shows  that  the  MFA  estimate  ranges  from  53?"  to 
78 % greater  than  present  data,  with  an  average  of  68 % greater. 
However,  the  MFC  estimate  ranges  from  1 % to  l6%  less  than  pre- 
sent data,  with  an  average  of  7-5'f3  less. 

Hall  and  Finger 

Fracture  measurements  of  flawed  holes  were  used  to  em- 
pirically relate  flaw  size  and  fracture  strength,  Kjc*  Their 
failure  criterion  assumed  the  form 


MFA  = (3-36)  (O.63)  = 2.12 


(10) 


MFC  =0.86  B( c/r ) 


(11) 
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KJc  = Ca  Vrrce  F(cg/r) 


(12) 


where  cg  is  an  effective  crack  length  corresponding  to  a thru- 
the-thickness,  Bowie  crack.  They  found  that  a value  for  C of 
0.87  correlated  well  with  their  data. 

For  present  estimates  the  form  of  Eq  12  was  used  to  esti- 
mate K.  Also  it  was  necessary  to  interpolate  the  curves  for 
ce  s^nce  Hall  and  Finger's  data  incorporated  a maximum  a/c  of 
one.  Using  C = 0.8?,  the  magnification  factor  expression  be- 
comes 


MF  = 0.87  F(ce/r)  (13) 

Table  I shows  good  agreement  with  present  data.  MFA 
estimates  range  from  11#  to  24#  below  present  data,  with  an 
average  of  16#  below.  MFC  estimates  range  from  5#  to  18#  be- 
low present  data,  with  an  average  of  11#  below. 

It  was  noted  that  a higher  value  of  C might  give  better 
agreement.  One  justification  for  a change  could  be  a signifi- 
cant difference  in  plastic  zone  sizes  of  the  testing  materials. 
Hall  and  Finger  used  two  test  materials,  2219-T87  aluminum  and 
5A1-2 . 5Sn(ELI)  titanium,  that  have  similar  plastic  zone  sizes, 
0.060  and  0.050  inch.  The  plastic  zone  size  of  PMMA  is  much 
smaller,  about  0.003  (Ref  18»33).  Accordingly,  a C value  of 
1.0  was  used  to  yield  modified  Hall  and  Finger  estimates. 

Table  I shows  very  good  agreement  with  present  data.  MFA 
estimates  range  from  3#  above  to  12#  below  present  data,  with 
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an  average  of  2.7^  below.  Excellent  agreement  is  shown  for  all 
values  of  a/T  except  the  higher  ones.  MFC  estimates » range  from 
2 fo  below  to  above  present  data,  with  an  average  of  2.5 °h  above. 
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VI.  Conclusions  and  Recommendations 


Conclusions 

1.  Fatigue  crack  growth  rate,  da/dn,  as  a function  of 
stress  intensity  factor,  K,  for  PMMA  under  constant  amplitude, 
cyclic  loading  can  be  characterized  by 

da/dn  = 0. 22138 ( 10 )"20K5* 591  (14) 

R-ratio  ( min/  max)  should  be  near  or  equal  to  zero.  For  this 
study  it  was  0.029. 

2.  Testing  at  1 or  2 cps  has  a negligible  effect  (<5%)  on 

results  in  the  range  of  crack  growth  from  0.15(10)""^  to  0.30(10)- 
inch/cycle . 

3-  Present  estimates  of  stress  intensity  factors  for 
corner  flaws  emanating  from  holes  tend  to  overestimate  K at 
the  hole  and  underestimate  K at  the  surface  for  flaw  shapes  in 
the  range  of  1.4.  A modified  Hall  and  Finger  empirical  rela- 
tion shows  good  agreement  with  present  data,  especially  for  K 
at  the  hole. 

4.  The  scatter  of  test  data,  particularly  at  early  growth 
stages,  is  probably  a result  of  the  quality  of  photographs  being 
very  sensitive  to  illumination  of  the  crack.  Scatter  in  the 
magnification  factor  plots  at  early  a/T  values  could  possibly 
be  a result  of  the  crack  not  reaching  its  natural  shape  by  that 
range  in  a/T.  A better  method  of  crack  initiation  would  probably 
yield  better  quality  pictures  and  more  efficient  examination  of 


early  growth. 

5.  The  fatigue  crack  growth  test  is  a useful  experimental 
technique  for  obtaining  stress  intensity  factor  calibrations, 
especially  for  complex,  three-dimensional  geometries. 

Recommendations 

1.  The  modified  Hall  and  Finger  relation  suggests  the 
possibility  of  using  the  fatigue  crack  growth  test  to  find  an 
empirical  relation  for  complex  flaw  geometries.  Further  crack 
growth  rate  tests  should  include  a study  to  find  a relation  be- 
tween stress  intensity  factor  and  crack  length  that  fits  the 
experimental  data  well  and  yields  conservative  estimates  of 
stress  intensities  that  can  be  used  either  for  approximation 

or  to  identify  trends. 

2.  Fatigue  crack  growth  tests  should  be  conducted  with 
various  specimen  thicknesses  and  hole  diameters  which  would 
yield  various  crack  shapes.  Thus  the  effect  of  crack  shape, 
a/c,  could  be  determined. 

3.  Tests  of  other  flaw  geometries  are  of  importance  and 
should  be  studied.  Geometries  such  as  the  completely  embedded 
flaw  at  a hole,  two  comer  flaws  at  a hole,  or  flaws  at  a hole 
with  partial  load  transfer  should  receive  primary  consideration. 
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Appendix 

Test  Data  and  Plots 
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Table  II 


Actual  Specimen  Dimensions 


Type  of  Test  Width, W Thickness, T Hole  Diam,2r 

(inches)  (inch)  (inch) 


Baseline 

7-95 

0.078 

0.7^2 

Baseline 

7. 98 

0.717 

0.740 

Comer  flaw 

7-98 

0.729 

0.747 

Corner  flaw 

7-97 

0.690 

0-753 

Corner  flaw 

7-98 

0.730 

0.750 

Corner  flaw 

7- 95 

0.720 

0.742 

Corner  flaw 

7-95 

0.698 

0.739 

Corner  flaw 

7-98 

0.704 

0.743 

Corner  flaw 

7-97 

0.704 

0.743 
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Fig.  15.  Magnification  Factor  at  the  Hole  for  a/c  = 1.2 
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Fig.  16.  Magnification  Factor  at  the  Surface  for  a/c  = 1.3 

' Q-trnr  ^ 


u~\ 

c'-n 


V 

o 

\ 

V 

CM 


.0 


.0 


00  ‘Z  09*1  02 ‘I  08 ’0  OV'O  00 '0° 

^013U J NOIldGIdINOdU 


44 


2 0*  24  0.36  0.48  0.60  0.72  0.84  0.96  ] . 08 

fl/T 

Fig.  17.  Magnification  Factor  at  the  Hole  for  a/c  = 1.3 
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Fig.  22.  Magnification  Factor  at  the  Surface  for  a/c 
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Fig*  23.  Magnification  Factor  at  the  Hole  for  a/c  = 1.6 
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Fig.  25-  Crack  Shape  vs.  a/T,  Test 
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Fig.  26  . Magnification  Factor  vs.  a/T,  Test  1 
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Fig.  2 7.  Stress  Intensity  Factor  vs.  c Crack  Length,  Test  1 
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Fig.  29  . Crack  Shape  vs.  a/T,  Test 
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Fig.  30.  Magnification  Factor  vs.  a/T,  Test  2 
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Fig.  31.  Stress  Intensity  Factor  vs.  c Crack  Length,  Test  2 
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Fig.  33.  Crac’K  Shape  vs.  a/T,  Test 
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Fig.  35-  Stress  Intensity  Factor  vs.  c Crack  Length,  Test  3 
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Fig.  36.  Crack  Length  vs.  Cycles,  Test  5 
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Fig.  37.  Crack  Shape  vs.  a/T,  Test 
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Fig.  39  • Stress  Intensity  Factor  vs.  c Crack  Length,  Test  5 
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Fig.  40.  Crack  Length  vs.  Cycles,  Test  6 
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Fig.  4l  . Crack  Shape  vs.  a/T,  Te 
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Fig.  43.  Stress  Intensity  Factor  vs.  c Crack  Length,  Test  6 
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Fig.  4 5.  Crack  Shape  vs.  a/T,  Test  8 
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Intensity  Factor  vs.  c Crack  Length,  Test  8 
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Fig.  ^9.  CracI:  Shape  vs.  a/T,  Tes 
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Fig.  51.  Stress  Intensity  Factor  vs.  c Crack  Length,  Test  9 
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. 447 

.506 

. 391 

.448 

23000 

.452 

.512 

.404 

.456 

24000 

• 4 63 

. 526 

.415 

.466 

25000 

.472 

.5  33 

. 426 

.477 

26000 

.502 

.562 

.451 

.505 
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28000 

29000 

31000 

33000 

35000 

37000 

38000 

39000 

41000 

43000 

4400  0 

45000 

46000 

4700  0 

48000 

50000 

51000 

52000 

53000 

54000 

55000 

56000 

57000 

58000 

5900  0 

60000 

61000 

62000 


AT  (INCH) 


AM(INCH)  A 8 ( INCH)  A A VG  ( INC  H ) 


. 499 

.564 

. 515 

.576 

.5  20 

.5  83 

. 5 42 

.609 

. 568 

.632 

.587 

. 653 

.621 

.687 

.622 

.688 

.628 

. 696 

.6  65 

.730 

. 718 

.777 

.714 

.779 

.733 

.802 

. 7 56 

.8  28 

.770 

.832 

. 791 

. 851 

.815 

. 881 

.833 

.890 

. 856 

.918 

.883 

.943 

.895 

.954 

.923 

.984 

.965 

1.012 

. 992 

1.  043 

1.022 

1.075 

1.058 

1.119 

1.  122 

1.170 

1.  168 

1.220 

1. 231 

1. 284 

.45  8 

.507 

.469 

.520 

. 431 

.529 

. 499 

.550 

.527 

.576 

.546 

.595 

. 580 

.629 

.585 

.631 

.601 

.642 

. 534 

.677 

♦ 67  b 

.724 

.63  3 

.725 

. ^Ob 

.7*7 

.734 

. 773 

. 7 37 

.780 

. 757 

.800 

. 787 

.827 

. 30C 

.841 

. 331 

.869 

. 854 

.893 

. 866 

.905 

. 90  0 

.936 

. 921 

• 96b 

.960 

.999 

. 99  0 

1.029 

1.  034 

1.070 

1.  097 

1.130 

1.  1*2 

1.177 

1.  193 

1.236 
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TEST  5 


NUMBER  OF  DATA  POINTS  26 


CYCLES  A < I 

1 

2000 
3000 
4000 
5000 
6000 
700  J 
8000 
9000 

10000  . 

11000 
1200  0 
13000 
14000 
15000 
16000 
17000 
18000 

18250  . 

18500 

18750 

1900  0 

19250 

19500 

19750 

19900 


C ( INCH) 

.128 
.131 
.134 
.141 
.148 
.154 
.161 
.167 
.175 
.184 
.198 
.208 
.220 
.234 
.253 
.272 
.291 
.315 
.322 
.326 
.335 
.348 
.359 
. 365 
.383 
.391 


NCH) 

135 

122 

125 

134 

143 

157 

166 

176 

190 

203 

216 

240 

255 

280 

310 

346 

390 

446 

465 

481 

496 

521 

538 

561 

592 

603 
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ST  6 

CLES 

1 

2000 

300  0 

<<000 

6000 

7000 

800  0 

9000 

11000 

12000 

1300  0 

14000 

15000 

1600  0 

17000 

1800  0 

19000 

20000 

20500 

21003 

21500 

22000 

22300 

22600 

22900 

23100 

23330 

23500 

23700 

23900 


NUMBER  OF  DATA  POINTS  30 


A (INCH) 

. 084 
, 097 
. 093 
. 107 
. 12  0 
. 142 
. 141 
. 154 
. 171 
. 184 
. 189 
. 204 
. 207 
. 237 
. 256 
. 282 
. 304 
, 346 
. 361 
. 376 
. 393 
. 427 
. 45  2 
. 46  3 
.480 
.496 
. 485 
. 502 
. 517 
.553 


C ( INCH ) 

.095 
.103 
.103 
.107 
.112 
.113 
.117 
.124 
.137 
.140 
.150 
.161 
.169 
.183 
. 192 
.202 
.218 
.238 
.248 
.256 
.270 
.290 
.302 
.310 
.322 
.331 
.343 
. 354 
.368 
. 383 
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TEST  7 


NUMBER  OF  DATA  POINTS  29 


CYCLES 

AM  (INCH) 

AB ( INCH) 

AA  VG  ( IN( 

1 

. 565 

. 432 

.499 

1000 

.550 

.428 

.489 

2000 

.578 

.452 

.515 

3000 

.601 

.481 

.541 

<*000 

.612 

. 493 

• 553 

5 C 0 0 

. 648 

.533 

.591 

6000 

. 642 

. 523 

. 583 

7000 

.649 

.536 

.593 

8 0 0 G 

. 693 

.579 

.636 

9000 

. 742 

.631 

.687 

10000 

. 745 

.641 

.69  3 

11003 

. 750 

.644 

.697 

1200  0 

. 773 

.674 

. 723 

13000 

. 777 

.682 

. 730 

14000 

.771 

.682 

.726 

15000 

.811 

.721 

. 766 

16000 

.816 

.728 

. 772 

17000 

.8  92 

.806 

. 849 

18000 

.918 

.832 

.875 

1900  0 

.912 

.832 

. 372 

20000 

. 950 

.862 

.906 

21000 

.952 

.864 

. 908 

22000 

.974 

. 898 

.936 

23000 

1.  014 

.938 

.976 

24000 

1.  043 

. 966 

1.  005 

26000 

1.  113 

1.038 

1.  075 

27000 

1.  163 

1.092 

1.127 

28000 

1.  183 

1.  115 

1.  149 

29000 

1.  262 

1.191 

1.  227 
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TEST  8 


NUMBER  OF  DATA  POINTS  34 


CYCLES  A(INCH)  C(INCH) 


1 

200  0 
3000 
4000 
500  0 
600  0 
700  0 
8000 
9000 
1000  0 
11000 
1 200  0 
13000 
14000 
15000 
16000 
17000 
18000 
19000 
20000 
21000 
22000 
22500 
23000 
23500 
24000 
24300 
24500 
24750 
25000 
25250 
25500 
25750 
26000 


. 118 
. 131 
. 133 
• 145 
. 146 
. 156 
. 168 
. 171 
.177 
. 183 
. 196 
.203 
. 212 
. 227 
. 248 
. 252 
. 274 
. 292 
. 312 
. 341 
. 362 
. 395 
.414 
• 440 
. 462 
. 495 
. 505 
.512 
.526 
.542 
.564 
. 580 
.600 
.612 


.132 
.130 
.129 
.133 
.129 
.127 
.129 
.131 
.132 
.139 
.144 
.152 
. 157 
.168 
.176 
.183 
.195 
.206 
.217 
.233 
.247 
.264 
.276 
.285 
.298 
.313 
.322 
.3  21 
.331 
.336 
.344 
.354 
.365 
.377 
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TEST  9 

NUM8ER  OF 

OATA  POINTS  45 

CYCLES 

A(INCH) 

C( INCH) 

16000 

.241 

.168 

16000 

• 251 

.177 

19000 

.261 

.184 

20000 

.279 

.192 

21000 

. 285 

.194 

21500 

.299 

• 202 

22000 

• 303 

.206 

22500 

. 307 

.210 

23000 

. 322 

.214 

23500 

. 331 

.223 

240  0 0 

. 342 

.229 

24500 

.353 

.231 

25000 

.359 

.239 

25500 

. 365 

• 240 

26000 

• 380 

• 246 

26500 

. 394 

• 256 

27000 

.409 

.263 

27300 

.409 

.267 

27600 

.424 

.274 

27900 

.434 

.279 

28200 

.441 

.284 

28500 

.452 

.292 

28800 

. 463 

.299 

29000 

.475 

• 303 

29200 

.471 

• 307 

29400 

.486 

• 314 

29600 

.494 

.319 

29800 

.502 

.323 

30000 

.510 

• 330 

30200 

.520 

• 338 

30400 

.529 

.340 

30600 

.539 

.347 

30800 

. 546 

• 354 

31000 

.572 

• 363 

31200 

.580 

.369 

31300 

• 586 

.371 

31400 

.599 

• 376 

31500 

• 606 

• 384 

31600 

.608 

• 387 
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CYCLES 


ft  (INCH) 


C (INCH) 


31700  .612  .392 
31800  . 620  » 394 
31900  .630  .401 
32000  .635  .406 
32050  *639  .411 
32100  .646  .412 


90 

☆ U.  S.  GOVERNMENT  PRINTING  OFFICE : 1975  — 657-630/58 


